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ABSTRACT 

Post-traumatic stress disorder (PTSD) is a serious mental health injury which can manifest after experiencing a 
traumatic life event. The disorder is characterized by symptoms of re-experiencing, avoidance, emotional numb- 
ing and hyper-arousal. Whilst its aetiology and resultant symptomology are better understood, relatively little is 
known about the underlying cortical pathophysiology, and in particular whether changes in functional connec- 
tivity may be linked to the disorder. Here, we used non-invasive neuroimaging with magnetoencephalography 
to examine functional connectivity in a resting-state protocol in the combat-related PTSD group (n = 23), and 
a military control group (n = 21). We identify atypical long-range hyperconnectivity in the high-gamma-band 
resting-state networks in a combat-related PTSD population compared to soldiers who underwent comparable 
environmental exposure but did not develop PTSD. Using graph analysis, we demonstrate that apparent network 
connectivity of relevant brain regions is associated with cognitive-behavioural outcomes. We also show that left 
hippocampal connectivity in the PTSD group correlates with scores on the well-established PTSD Checklist (PCL). 
These findings indicate that atypical synchronous neural interactions may underlie the psychological symptoms 
of PTSD, whilst also having utility as a potential biomarker to aid in the diagnosis and monitoring of the disorder. 
® 2014 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY-NC-ND license 

(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Post-traumatic stress disorder (PTSD) is a mental health problem, 
characterized by anxious and depressive features, which develops 
after exposure to a traumatic life event. It is placed in the Diagnostic 
and Statistical Manual of Mental Disorders (DSM-V) trauma and 
stressor-related disorders category, and PTSD is principally comprised 
of four symptom clusters: re-experiencing; avoidance; emotional 
numbing; and hyperarousal (American Psychiatric Association, 2013). 
Incidence of the disorder in the general population is around 5-10% 
(Kessler et al., 2005), and is thought to be much higher in some military 
populations returning from recent combat deployments in hazardous 
regions (Richardson et al., 2010). Knowledge remains scant, however, 
regarding the neurobiological basis of PTSD, limiting our understanding 
of the disorder and hampering the search for reliable biomarkers. 
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StiTictural and functional neuroimaging studies of PTSD using mag- 
netic resonance imaging (MRl and fMRl) and positron emission tomog- 
raphy (PET) have reported atypical neuroanatomy and differential 
activation patterns in a number of cortical and subcortical structures 
(Hull, 2002). Known to play a role in episodic memory, the hippocampi 
show both structural and functional abnormalities (Bremner et al., 
2003 ) . Additionally, the amygdalae ( Etkin and Wager, 2007 ) , ventrome- 
dial prefrontal (Gold et al., 2011), and the dorsal anterior cingulate 
cortices (Shin et al., 201 1 ) have been reported to show atypical function 
in the disorder (for an extensive biological review of PTSD, see Pitman 
et al.,2012). 

In recent years, advances in the analyses of resting-state networks 
and the functional connectivity among brain regions have allowed re- 
searchers to map ongoing and spontaneous communication required 
for the temporal coordination of cognitive and sensory processing 
(Damoiseaux et al., 2006). In particular, changes in functional connec- 
tivity defined by ongoing oscillatory synchrony (Wang, 2010) have 
proven to be useful in mapping cortical pathophysiology, thought to un- 
derlie a number of neurophysiological disorders (Tewarie et al., 2013), 
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including magnetoencephalograpliic (MEG) studies of PTSD. Tlie tecli- 
nique of MEG allows spatiotemporal patterns of brain activity to be 
mapped with millisecond precision, elucidating functional brain chang- 
es on time scales to which fMRl is blind (Hari and Salmelin, 2012); fMRl 
only images ultra-low frequencies. Furthermore, MEG benefits from 
being able to directly record ongoing brain activity, as opposed to 
changes in the haemodynamics associated with neural firing. 

MEG studies of PTSD suggest enhanced slow wave generators in left 
temporal regions, and decreased oscillations in parieto-occipital cortex 
are related to PTSD (Kolassa et al., 2007). Georgopoulos et al. (2010) 
suggested that patterns of abnormal synchronous oscillations could dif- 
ferentiate PTSD from control subjects, particularly poor communication 
between right temporo-parietal areas and other brain regions (Engdahl 
et al., 2010). This group also linked decorrelations in small networks, 
most evident in the right superior temporal gyrus, with resilience to life- 
time trauma in control veterans, but not those with PTSD (James et al., 
2013). The above studies collectively suggest that abnormal coherent 
brain oscillations might contribute to symptoms of the disorder. 

Prior studies, however, have not investigated frequency-specific in- 
teractions in source-resolved networks, or their association with expo- 
sure to stressful stimuli or symptom severity. Here we used source- 
analysed MEG and graph theoretical analysis to test the hypotheses 
that veterans with PTSD would express atypical resting-state network 
synchrony; that these atypical inter-regional interactions would be 
exacerbated by exposure to stressful, combat-related imagery; and 
that an aberrant organization of neurophysiological networks is associ- 
ated with the severity of PTSD symptoms and associated cognitive- 
behavioural outcomes. We selected strength and degree as our graph 
properties of interest, as they most directly correspond to network 
hyperconnectivity or hypoconnectivity when comparing two popula- 
tions, and thus correspond most closely with our hypothesis that vet- 
erans with PTSD would show alterations in the intensity of network- 
level neurophysiological interactions. Specifically, node strength repre- 
sents the weighted magnitude of the relation between a seed node and 
the network, which is the sum of all connections, and node degree indi- 
cates the number of connections between a given node and other nodes 
in the network (Bullmore and Sporns, 2009), which would capture dif- 
ferential network synchrony between veterans with PTSD and their 
matched controls. 

2. Methods and materials 

2.1. Participants 

MEG data were recorded from 23 Canadian Armed Forces soldiers, 
who deployed in support of the Afghan mission and were subsequendy 
diagnosed with PTSD (all male, mean age = 37.4, SD = 6.8, age range 
22-48). Twenty-one soldiers (all male, mean age = 33.05, SD = 5.26, 
age range 18-45) who also participated in the Afghan mission but did 
not develop PTSD were recruited as a control group. 

Participants were included in the PTSD group if they met the follow- 
ing criteria: they have a diagnosis of combat-related PTSD from an oper- 
ational trauma stress support centre (OTSSC); PTSD symptoms were 
present from 1 to 4 years prior to participation in the study; they were 
engaged in regular mental health follow-up; and they had moderate 
or greater severity on the PTSD Checklist (PCL > 50). The diagnosis 
was determined by a psychiatrist or psychologist specializing in 
trauma-related mental health injuries and conducted through a com- 
prehensive, semi-structured intei"view based upon DSM-IV-TR diagnos- 
tic criteria (American Psychiatric Association, 2000), along with 
Canadian Armed Forces (CAF) standardized psychometric testing. 
Interview-based clinician diagnosis of mental disorders is considered 
superior to pen and paper or self-administered screening methods. All 
participants in the PTSD group were recruited from one of the CAF 
OTSSCs, which are centres of excellence for the diagnosis and treatment 
of trauma-related mental health injuries. There were usually more than 



one DSM-IV-TR 'Al' stressor-related criteria (American Psychiatric As- 
sociation, 2000) identified as a traumatic event contributing to the de- 
velopment of PTSD (direct personal experience of an event that 
involves actual or threatened death or injuiy), with diagnosis related 
to operational exposure. Control soldiers were matched on rank, educa- 
tion level, handedness and militaiy experience. 

Additional inclusion criteria applied to both groups included: no his- 
tory of a traumatic brain injury (TBI), screened by a psychiatrist through 
a review of their electronic health record, telephone inte™ew, and ad- 
ministration of the Defence and Veteran's Brain Injuiy Centre (DVBIC) 3 
item screening tool; English-speaking and able to understand task 
instructions and give informed consent. Exclusion criteria included fer- 
rous metal inside the body that might be classified as MRl contraindica- 
tions or items that might interfere with MEG data acquisition; presence 
of implanted medical devices; seizures or other neurological disorders, 
or active substance abuse; certain ongoing medications (anticon- 
vulsants, benzodiazepines, and/or GABA antagonists) known to directly 
or significantly influence electroencephalographic (EEG) findings. As 
this was a naturalistic sample, however, all PTSD patients were on 
evidenced-based psychotropic medication(s), such as selective seroto- 
nin reuptake inhibitors (SSRls), serotonin-norepinephrine reuptake 
inhibitors (SNRls) and Prazosin. 

All participants underwent cognitive-behavioural testing in addition 
to the MEG resting-state scan. These assessments included: estimates of 
IQ from the Wechsler Abbreviated Scale of Intelligence (WASl); the 
Alcohol Use Disorders Identification Test (AUDIT); Conner's Attention- 
Deficit Hyperactivity Disorder Test; the Generalized Anxiety Disorder 
7 (GAD-7) test; Patient Health Questionnaire (PHQ9); and the Post 
Traumatic Stress Disorder Checklist (PCL). Within the PTSD group, 
there were significant rates of co-morbid mental disorders, including 
major depressive disorder (MDD; 74% of the PTSD group). These find- 
ings are consistent with prevalence rates established through large 
scale studies in military populations (Gather et al., 2012). 

2.2. Procedure and MEG data acquisition 

Resting-state MEG data were collected in two separate runs. Partic- 
ipants were supine and instructed to rest with eyes open and maintain 
visual fixation on an x within a circle on the screen. Following the first 
resting-state run, the participants completed a number of imaging 
protocols, with a memory-related paradigm containing trigger images, 
such as scenes of traumatic events (e.g. battlefield casualties) intermixed 
with neutral images, an emotional faces task, and a verbal task that 
contained neutral as well as salient trigger words (such as 'Kandahar' 
and 'grenade'). They then completed a second resting-state run. We 
expected the affective stimuli to induce arousal and attentional mecha- 
nisms, and perhaps differentially activate the PTSD group compared 
to the controls. We refer to the initial scan as the pre-triggering 
resting-state run and the subsequent acquisition as the post-triggering 
resting-state run. 

MEG data were collected inside a magnetically-shielded room on 
a CTF Omega 151 channel system (CTF Systems, Inc., Coquitlam, 
Canada) at The Hospital for Sick Children, at 600 Hz for 300 s per 
resting-state mn. Throughout the run, head position was continuously 
recorded by three fiducial coils placed on the nasion, and left and right 
pre-auricular points. 

After the MEG session, anatomical MRl images were acquired using 
the 3 T MRl Research scanner ( Magnetom Tim Trio, Siemens AG, Eriangen, 
Germany) in a suite adjacent to the MEG. Structural data were obtained as 
Tl-weighted magnetic resonance images using resolution 3D MPRAGE 
sequences (repetition time [TR] = 2300 ms; echo time [TE] = 2.9 ms; 
flip angle [FA] = 9°; field-of-view [FOV] = 28.8 x 19.2 cm; 256 x 256 
matrix; 192 slices; 1 mm isovoxel) on a 12-channel head coil. MEG data 
were coregistered to the MRl structural images using the reference 
fiducial coil placements. A multi-sphere head model was constructed for 
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each individual and ttieir brain space was normalized to a standard 
Montreal Neurological Institute (MNI) brain using SPM2. 

2.3. MEG data processing 

2.3.1. Seed definition and virtual electrode recording 

MEG data were band-pass filtered offline at 1-150 Hz, a notch filter 
applied at 60 Hz (8 Hz bandwidth) and a third-order spatial gradient 
environmental noise-cancellation applied to the recording. A priori 
sources (seeds) of interest in cortex and sub-cortical regions were iden- 
tified from the Automated Anatomical Labeling (AAL; Tzourio-Mazoyer 
et al., 2002) atias giving 90 locations for time-series to be extracted and 
analysed. Broadband time-series (Virtual electrodes') from these voxels 
were reconstructed using a vector beamformer on the basis of the 90 
AAL coordinates for each subject and filtered into five classical EEG 
bandwidths for further analysis: Theta (4-7 Hz), Alpha (8-14 Hz), 
Beta (15-30 Hz), 'low' Gamma (30-80 Hz), and 'high' Gamma 
(80-150 Hz). 

Beamformers are a type of spatial filter used to suppress signals from 
unwanted noise sources, whilst being optimally sensitive to activity in a 
given brain location (in this particular case, the AAL seed locations). 
Individual weight vectors are applied to each sensor measurement and 
summated to give estimated source activity to a particular cortical seed 
location (Quraan and Cheyne, 2010). Additionally, MEG beamformers 
are effective at suppressing ocular artefacts generated by eye movements, 
which are a particular problematic source noise in EEG, and non-ocular 
artefacts, such as cardiac and muscle activity (Muthukumaraswamy, 
2013). Thus, beamforming removes artefacts, not requiring the rejection 
of data. 

2.3.2. Assessing fiinctional connectivity: weighted phase lag index 

Each of the 5 band-pass filtered waveforms were then submitted to a 
functional connectivity analysis, using the weighted phase lag index 
(wPLl; Vinck et al., 2011). The instantaneous phase of each sample 
from the filtered time-series was calculated using the Hilbert transform. 
The degree of phase synchronization between all pairwise combinations 
of the seeds was computed using the wPLl, which is based on the mag- 
nitude of the imaginary component of the cross-spectrum (Lau et al., 
2012). Ranging between 0 and 1, these values quantify the phase syn- 
chrony between two cortical/sub-cortical sources, referred to as func- 
tional connectivity. 

2.3.3. Statistical analysis 

Adjacency matrices with wPLI values acting as edge weights for all 
sources were constructed, which resulted in a 90 x 90 [x5 frequencies] 
matrix of weighted undirected graphs for each participant. Adjacency 
matrices were then divided into groups and inferential statistics inves- 
tigating group differences for edge weights were implemented using 
non-parametric permutation testing, which do not require that the 
data distributions be normal, using the Network Based Statistic (NBS) 
Toolbox (Zalesl<y et al., 2010). False positives due to multiple compari- 
sons were controlled for using the false discovery rate (FDR; corrected, 
p < 0.05), across the whole 90 x 90 matrix. Network measures of node 
strength were used to assess the importance of a node within the net- 
works and these measures were obtained using Brain Connectivity 
Toolbox (Rubinov and Sporns, 2010). Brain networks were visualized 
using BrainNet Viewer (Xia et al., 2013). Further behavioural corre- 
lation analyses were conducted using MATLAB Statistics Toolbox 
(The Mathworks, Inc.). 

3. Results 

Measures for the cognitive-behavioural test batteries and the 
characteristics of the soldiers with and without PTSD are presented 
in Table 1. 



Table 1 

Cognitive-behavioural assessment measures for PTSD and conti'ol participants. 



PTSD Control Test statistic 



n 


23 


21 




WASl 


108 (14.09) 


116.81 (13.67) 


t = 2.10, d/= 42, p = 0.0420 


AUDIT 


8.70 (7.67) 


5.10 (3.91) 


t= -1.40, (i/=42. 








p = 0.170 


Conner's 


24.09 (9.57) 


7.19 (4.99) 


t= -7.24, d/= 42 p< 0.001 


CAD-7 


15.04 (4.41) 


1.95 (1.99) 


t= -12.50, d/= 42, 








p< 0.001 


PHa9 


16.87 (5.09) 


1.86 (2.31) 


t= -12.39, d/= 42, 








p < 0.001 


PCL 


63.95 (7.81) 


NA 





WASl, Wechsler Abbreviated Scale of Intelligence; AUDIT. Alcohol Use Disorders Identifi- 
cation Test; Conner's, Attention-Deficit Hyperactivity Disorder Test; GAD-7, Generalized 
Anxiety Disorder 7; PHQ9, Patient Health Questionnaire; PCL, Post Traumatic Stress Disorder 
Checklist 



3. 3 . Functional connectivity analyses: pre- and post-triggering resting-state 
networks 

Increased high gamma (80-150 Hz) resting oscillatoiy network syn- 
chrony was observed in the PTSD soldiers compared to control soldiers 
in the pre-triggering resting-state run (p < 0.05, FDR-corrected; Fig. 1 A). 
This network comprised 34 nodes of varying degree and strength, local- 
ized predominantly in the left hemisphere, with connectivity differ- 
ences in important nodes noted in the left frontal, left temporal, and 
right parietal regions. For the other frequency bands, the analysis 
revealed only the occasional significant edge connecting two nodes 
(p < 0.05), but these patterns are inconclusive in terms of revealing 
any large-scale network connectivity differences, and were not investi- 
gated further. 

All seed regions found to have significantly different connectivity 
patterns as a result of the permutation analysis (FDR-corrected, 
p < 0.05) are listed in Fig. 2A, first ordered by hemisphere, and then by 
lobe. Fig. 2A illustrates the magnitude of the difference between the 
groups in node degree (in other words, the number of significant con- 
nections or edges between nodes). In Fig. 2B individual group-mean 
measures of node strength are shown (the sum of edge weights for a 
particular node; PTSD soldiers denoted by the black bars, control 
soldiers by the white bars). Important, high-strength nodes with a 
degree difference greater than 3 are highlighted by the asterisk in 
Fig. 2A. None of the other frequency bands displayed significant atypical 
network connectivity in the between-groups pre-triggering resting- 
state analysis. 

A within-group comparison of pre- versus post-triggering resting- 
state networks revealed a similar pattern of network connectivity 
changes for the control soldiers after viewing trigger images. In the 
post-resting-state they displayed enhanced high gamma band connec- 
tivity compared to the pre-triggering resting-state run, with increases 
in connectivity confined mostly to the left hemisphere (p < 0.01, FDR- 
corrected; Fig. IB, bottom row); the most highly connected nodes 
were the left precuneus (node degree = 4) and left amygdala (node 
degree = 4). The PTSD soldiers showed relatively minor increases 
in high gamma network connectivity between the pre- and post- 
triggering resting-state sessions (p < 0.01, FDR-corrected; Fig. IB, top 
row), with enhanced inter-hemispheric synchrony involving the left 
middle frontal gyrus (node degree = 2) and the right supramarginal 
gyrus (node degree = 2). 

Finally, a between-group comparison in the post-triggering resting- 
state run revealed only minor network differences after viewing the 
trigger images, with the right orbital-frontal cortex and right amygdala 
showing enhanced connectivity in the PTSD group compared to the 
controls (p < 0.05, FDR-corrected; Fig. IC). This suggests viewing 
emotionally-salient stimuli brings about network connectivity changes 
to the control group that resemble the atypical resting-state synchrony 
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Pre-triggering Resting-State Networks — Between Groups 



A 




FDR-corrected, p < 0.05 



Post- versus Pre-triggering Resting-State Networl<s — Within Groups 




FDR-corrected, p < 0.01 



Post-triggering Resting-State Networl<s — Between Groups 




FDR-corrected, p < 0.05 
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Fig. 2. Summary of brain regions identified as being hyperconnected in the PTSD group compared to controls in the pre-triggering resting-state run derived from the NBS analyses (p < 0.05, 
FDR-corrected). (A) The list of brain regions that were found to be significant is initially divided by hemisphere, then by lobe, and then finally by degree in ascending order. Nodes marked 
with an asterisk (*) indicate high-strength, critical network nodes (degree > 4) that were used in the cognitive-behavioural correlation analyses. (B) A between-group comparison of the 
mean within-network node strengths (sum of the wPLl) values. 



shown in the PTSD group, such that these otherwise-similar popula- 
tions who underwent comparable experiential events show broadly 
comparable network connectivity patterns after the triggering stimuli, 
with a small increase in connectivity between the right amygdala and 
right orbitofrontal gyrus. 

32. Seed connectivity correlates of cognitive-behavioural outcomes 

To test the hypothesis that aberrant resting-state synchrony is associ- 
ated with cognitive-behavioural outcome and symptom severity, correla- 
tions of node strength against test measures were performed (Fig. 3). To 
understand inter-subject variability in the neuropsychological data, both 
groups were combined as the participants were experientially-similar in- 
dividuals, likely to exhibit cognitive-behavioural and affective tendencies 
along a spectnim. 

Network strength values were not normally distributed and there- 
fore correlation analyses were performed using Spearman's rank order 
correlation. These analyses were restricted to high degree nodes that 
were found to be significant following the analysis, which provided 6 
seeds/region of interest (ROI) in 3 lobes to explore: left frontal mid- 
orbital cortex in frontal regions; the left middle temporal gyrus, left pos- 
terior cingulum, left hippocampus, and left inferior temporal gyrus in 



the temporal lobe; and finally the right precuneus in parietal cortex 
(* in Fig. 2A denotes the nodes used in the correlation analysis). Impor- 
tantly, these nodes were found in regions that are often implicated in 
the PTSD literature (Pitman et al., 2012). 

Measures of anxiety (based on the GAD-7 score) and depression 
(from the PHQ9) were found to significantly correlate (p > 0.05, FDR- 
corrected) with the atypical network gamma connectivity strength in 
the left inferior temporal (Fig. 3 A) and left orbitofrontal regions 
(Fig. 3B). Furthermore, correlations for depression scores were also 
found to correlate with connectivity in the left hippocampus (Fig. 3C) 
and right precuneus (Fig. 3D). Finally, scores from the PTSD Checklist 
(PCL; characterizes PTSD symptomology) were found to correlate high- 
ly with gamma network strength in the left hippocampus in the PTSD 
Soldiers (Fig. 3C). 

4. Discussion 

This study provides the first evidence of frequency-specific alter- 
ations in source-resolved network synchrony in PTSD. In particular, 
combat-related PTSD soldiers display long-range hyperconnectivity in 
the high-frequency gamma band compared to matched controls who 
did not develop PTSD. We uniquely demonstrate that exposure to 



Fig. 1. Network connectivity maps in the high gamma band, and inter- and intra-regional adjacency matrices. (A) Contrasts of pre-triggering resting-state networks (RSN) for the PTSD 
versus control soldier participants in the 80-150 Hz band, showing hyper connectivity in the PTSD group compared to controls (p < 0.05, FDR-corrected) with node size scaled to degree. 
Mean wPLl region x region adjacency matrbc for PTSD minus control group, with each element of the graph showing vertices depicting hemisphere-divided lobes, such that the top left and 
bottom right quadrants denote intra-hemisphere connectivity (left and right, respectively); the top right and bottom left quadrants show inter-hemispheric lobe connectivity. Regions 
preceded by L denote left hemisphere, R denotes right hemisphere; F, frontal; T. temporal, SC, sub-cortical; P, parietal; 0, occipital. (B)Top row, within-group pre- versus post-triggering 
RSN for the PTSD participants in the 80-1 50 Hz band showing increased connectivity in the post-triggering RSN (p < 0.01, FDR-corrected). Bottom row, same as above except for the control 
soldier participants. (C) Post-triggering RSN for PTSD versus control groups (p < 0.05, FDR-corrected). 
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Fig. 3. Networif connectivity identified in pie-tiiggering RSN for FTSD versus contiol groups, and within-networif node connectivity (wPLl; undirected weiglited stiengtli measures) cor- 
related with cognitive-behavioural assessment outcomes (p < 0.05, FDR-corrected). (A) IMixed/combined ITSD and control group left inferior temporal cortex node strength correlated 
against GAD-7 (top plot, r (42) = 0.48, p = 0.002), and PHQ9 {bottom plot r (42) = 0.53, p = 0.001 ). (B) Left frontal orbital strength versus GAD-7 (top plot r (42) = 0.36, p = 
0.032) and PHQ9 {bottom plot r (42) = 0.42, p = 0.022). (C) In the bounded box, left hippocampus PTSD Group-only connectivity strength versus PCI {top plot r (18) = 0.58, p = 
0.042). Below, node sti'ength versus PHQ9 {bottom plot r (42) = 0.49, p = 0.001 ). (D) Right precuneus strength versus PHQ9 (r (42) = 0.56, p < 0.001 ). 



Stressful combat-related imagery can exacerbate atypical network con- 
nectivity associated with PTSD in military veterans. Furthermore, this 
study is the first to provide evidence that altered network topologies 
of critical brain regions are associated with cognitive-behavioural 
outcomes and symptom severity in PTSD, suggesting frequency- and 
region-specific alterations in neurophysiological synchrony may serve 
as a biomarker for PTSD. 

In contrast to an experientially-similar control soldier population, 
the PTSD group show significant network connectivity abnormalities 
in regions often implicated in PTSD, such as the episodic memory- 
related hippocampal (Bremner et al., 2003; Etkin and Wager, 2007; 
Pitman et al., 2012) and frontal regions (Gold et al., 2011; Shin et al., 
2004). Additionally, we found connectivity differences in temporal 
(Kolassa et al., 2007; Engdahl et al., 2010; James et al., 2013) and parietal 
regions (Kolassa et al., 2007) which are largely consistent with those re- 
ported in the literature (Pitman et al., 2012). 

One striking result of this study shows that following tasks involving 
emotionally-salient trigger images and words, the triggering brings the 
brain oscillations of veterans, who never developed PTSD, in line with 
those with the disorder. In other words, the difference between pre- 
and post-triggering resting-states in the control soldier group is similar 
to the connectivity differences between the groups in the pre-triggering 
resting-state. This suggests that trauma events affect connectivity simi- 
larly in both groups, but that individuals with PTSD are perhaps unable 
to return to a baseline state of ongoing, inter-regional oscillatory 
synchrony. 



We speculate that this hyperconnectivity involving the left 
hippocampus, temporal and frontal regions reflects some of the primary 
positive symptoms of PTSD, principally comprised of disturbing mental 
imagery and chronic hyper-arousal, possibly caused by re-experiencing 
and re-imagining of traumatic events, as well as heightened vigilance to 
aversive stimuli. This view is supported by findings from human intra- 
cranial recordings indicating that hippocampal-cortical gamma syn- 
chronization is associated with the formation of episodic memories 
(Fell et al., 2001 ) and in states of vigilance (Llinas and Steriade, 2006), 
both of which are heavily implicated in PTSD symptomology. This inter- 
pretation is further supported by an fMRl study that showed hippocam- 
pal activation in a fear response study correlated with re-experiencing 
and hyperarousal symptom cluster scores (Sripada et al., 2013) on a 
subscale from the Clinician-Administered PTSD Scale (CAPS). Moreover, 
fMRl resting-state studies of combat veterans have reported increased 
connectivity between the amygdala and insula in those with PTSD 
(Rabinak et al., 201 1 ). These changes were interpreted to be related to 
threat perception and modulated fear responses, cognitive states that 
are thought to be atypical in PTSD (Rabinak et al., 201 1 ). Similarly, we 
also demonstrated increases in connectivity between the right amygda- 
la and orbitofrontal cortex in the post-triggering resting-state results, 
which we believe are associated with an abnormal conditioned fear re- 
sponse ubiquitous in PTSD (Shin et al., 2004), which is also consistent 
with prior MEG reports of increased gamma event-related synchroniza- 
tion of the right amygdala during perception of emotional stimuli (Luo 
et al., 2009). We consider that this particular obsereation may be related 
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to the chronic heightened state of arousal that is symptomatic of the 
disorder. 

Despite similarities between our findings and the results of prior 
fMRl investigations, there are also discrepancies between the current 
study and other fMRl resting-state results in PTSD. A number of groups 
report simultaneous increases and decreases in fMRl functional connec- 
tivity in distinct networks in PTSD patients (Yin et al., 2011; Jin et al., 
2013). Yin et al. (201 1 ) found increased connectivity between the thal- 
amus and bilateral inferior frontal and left middle frontal gyri, as well as 
the left inferior parietal and right superior parietal regions. These re- 
gions have also been shown to be preferentially activated in PTSD pa- 
tients during a memoiy encoding task. The left middle frontal gyrus, 
in particular, (Bremner et al., 2003), is an area implicated in memory 
recollection. As well, the inferior parietal lobule, has been shown to be 
differentially active during presentation of emotionally-salient cues 
(Pagani et al., 2010). Jin and colleagues report atypical connectivity 
between the right PCC and insula, which they interpreted as underlying 
altered coupling between memory and perception, subserving the re- 
experiencing of traumatic episodic events (Jin et al., 2013). Both groups 
reached similar conclusions regarding the nature of these atypical 
network interactions, speculating that this enhanced connectivity may 
be related to positive symptom components of PTSD (such as chronic 
hyperarousal and traumatic memories), which corroborate the inter- 
pretation of our results. 

Yin and colleagues, however, also reported decreased connectivity be- 
tween the thalamus and right medial frontal gyius and left ACC, as well as 
a negative relation between CAPS scores and connectivity between the 
thalamus and right precuneus (Yin et al., 201 1 ). Similarly, Jin et al. report- 
ed decreased connectivity between the right amygdala and left middle 
frontal gyrus in PTSD which was negatively associated with CAPS scores 
(Jin et al., 2013); other studies also report reduced connectivity in 
resting-state networks, particularly the default mode network, as well 
as in the amygdala and hippocampal gyri (Bluhm et al., 2009). Taken to- 
gether, these effects are more difficult to reconcile with the results of the 
present study, although it should be noted that these relations existed in 
different neural systems, measured using distinct graph theoretical met- 
rics. Moreover, these previous resting-state results were obtained using 
different imaging modalities with dissimilar mechanisms of signal gener- 
ation with differing spatial and temporal resolutions (see Hall et al., 2013 
for a review of the relationship between MEG and fMRl). 

One important point to note is that our ability to interpret similari- 
ties and differences between the results of the present studies and 
prior findings with fMRl relies on our emerging understanding of the 
neurophysiological basis of the BOLD signal. It is well known that neuro- 
nal firing rate, high-frequency LFP power changes and the BOLD 
response are tightly coupled (Mukamel et al., 2005; Nir et al., 2007). 
Nonetheless, the relation between functional connectivity defined by 
BOLD amplitude covariations and high-frequency oscillatory dynamics 
are less well understood. Some report that BOLD amplitude covariations 
are positively correlated with gamma amplitude envelope changes 
(Tagliazucchi et al., 2012) and that ongoing BOLD fluctuations represent 
changes in cortical excitability and neural firing (Hiltunen et al., 2014). 
Conversely, others suggest only low-frequency oscillatory components 
influence BOLD cross-correlations (Wang et al., 2012), with a negligible 
contribution from ongoing gamma activity. 

Some researchers have highlighted the potential mechanistic and 
functional differences between 'intrinsic coupling modes' (ICMs; Engel 
et al., 2013) — two distinct modes of operation in oscillatory brain net- 
works. The first, 'envelope ICMs', arise from fluctuations in the ampli- 
tude envelopes (or power) of disparate neuronal populations, and 
their signal cross-correlations which reveal the degree to which they 
are communicating. A second, distinct type of ICM is the 'phase ICM', 
which is defined by the phase synchronization between the band- 
limited oscillatory signals emanating from different neural populations, 
in which communication is purported to be enabled through coherence 
(Fries, 2012). Here, we examined phase ICMs. 



Despite a lack of research examining the interplay between these 
two modes of action, it nevertheless has been proposed that these 
network-mediating mechanisms operate largely independently of one 
another (Engel et al., 2013). Moreover, changes in the degree to which 
coherently oscillating, or phase synchronized, neuronal groups couple 
might not necessarily incur increased metabolic demands, which 
would otherwise result in increases in both the ongoing ultra-slow 
BOLD fluctuations and/or amplitude envelope covariations seen in 
other studies. Accordingly, although experimental evidence indicates 
that BOLD correlations may bear some correspondence to inter-regional 
covariance of envelope fluctuations, the relation between BOLD correla- 
tions and inter-regional phase synchrony, such as that measured in the 
present study, remains more poorly understood. 

Finally, we show a significant correlation between left hippocampal 
node strength within the aberrant functional network and PTSD symp- 
tom severity on the PTSD Checklist in the soldiers with PTSD. This novel 
finding significantly extends previous research on ongoing local neural 
synchrony/cortical oscillations (Kolassa et al., 2007) by exploring the 
network interactions with behavioural function. We have gone beyond 
studies that have shown sensor-space (non-source resolved) measures 
of neural interactions in right temporal regions might be linked to PTSD 
and measures of resilience Qames et al, 2013), and shown that source- 
resolved activity in an important episodic-memory related structure 
correlates significantly with symptom severity in PTSD. 

4.3. Conclusions 

This study provides the first evidence that militaiy veterans who 
underwent similar environmental exposure but are differentiated by a 
post-deployment PTSD diagnosis, show differences in regional connectiv- 
ity in high-frequency oscillatory synchrony. Additionally, we uniquely 
demonstrate that 'high-level' areas associated with executive function 
and memoiy (the frontal and temporal lobes in particular), and the hip- 
pocampus, which is associated with episodic memory, show significant 
cortelations between high gamma neural networks and behavioural out- 
comes and symptom severity in PTSD. 
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